/*(PARI) _R. J. Cano_, Mar 6 2013*/
fsf(e,x)=((0<=e)&&(e<=4))*if(e==0,x,1)*if(e==1,(x*(x+1))\2,1)*if(e==2,(x*(x+1)*(2*x+1))\6,1)*if(e==3,(x*(x+1))^2\4,1)*if(e==4,(x*(x+1)*(2*x+1)*(3*x^2+3*x-1))\30,1); /*Fast summation facility*/
sievesum(n,k,e)=fsf(e,n)-2*(k^e)*fsf(e,(n\k));
A181482(n)=sievesum(n,3,1);
A213203(n)=sievesum(n!,n,1);